	PRO TVS_VEC,dir,time,Nproc
;TVS of J_z in the background (white is positive and black is negative). 
;Also shows the proton flow and the electron flow. 
;The largest arrow lengths correspond to the Alfv\`en speed.

Nlev=10 & Nreb=32
leer_input,dir,N,nu,eta,eps,epse,dt,nstep,nstg,nstdat,nstk,nstp
leer,dir,time,Nproc,a,f,b,u,jj,ww
print, eps
x=2.*!pi*dindgen(N)/N & y=x
dx=x(1)
fe=f-eps*b
comcampos,a,bx,by,x,y ; Compute magnetic field components
comcampos,f,ux,uy,x,y ; Compute velocity field components
comcampos,fe,uex,uey,x,y ; Compute electron flow components
loadct,6
;set_plot,'ps'
;device,file='fig_proton.eps',xs=12,ys=12,/color,/encapsulated
window,0,xs=N,ys=N,tit='Proton'
tvscl,jj
velovect,rebin(ux,Nreb,Nreb),rebin(uy,Nreb,Nreb),rebin(x,Nreb),rebin(y,Nreb),col=0,len=3*max(ux),chars=.01,/noerase
;device,/close
;device,file='fig_electron.eps',xs=12,ys=12
window,1,xs=N,ys=N,tit='Electron'
tvscl,jj
velovect,rebin(uex,Nreb,Nreb),rebin(uey,Nreb,Nreb),rebin(x,Nreb),rebin(y,Nreb),col=0,len=3*max(uex),chars=.01,/noerase
;device,/close
;set_plot,'x'
end
